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(57) Abstract: The invention relates to a 
correct inversion formula for tomographic 
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of the object (E), i.e. including translations, 
rotations and homotheties that are comparable 
to uniform dilatations. The invention also 
relates to generalisations relating to other 
deformation situations of the object studied 
and other radiations. 

(57) Abrege : Une formule d'inversion cor- 
recte de projections de tomographic est pro- 
posee en supposant une deformation affine de 
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tions, des rotations et des homotheties compa- 
rables a des dilatations uniformes. Des gene- 
ralisations a d'autres situations de deformation 
de l'objet 6tudie et d'autres rayonnements sont 
possibles. 
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PROCEDE DE RECONSTRUCTION D 1 UNE IMAGE 
TOMOGRAPHIQUE PAR UNE METHODE ANALYTIQUE COMPRENANT UNE 
MODELI SATION AMELIOREE DU MOUVEMENT DE L 1 OB JET . 

5 DESCRIPTION 

Le sujet de cette invention est un proc6d6 de 
reconstruction d'une image tomographique, du genre 
analytique et dans lequel on recourt a une modelisation 

10 perfectionnee du mouvement de l'objet afin de r6duire 
ies artefacts de 1 T image. 

La reconstruction d f images en tomographie 
implique 1 'utilisation d'un rayonnement traversant 
l'objet. Dans de nombreux cas, le rayonnement est un 

15 rayonnement d f irradiation partiellement attenue par les 
points de l'objet qu'il traverse ; il arrive aussi que 
les lignes de rayonnement soient fictives et 
correspondent simplement a des lignes de collimation 
des d6tecteurs, qui enregistrent une emission de 

20 particules dans l'objet. Les premiers de ces proc§des 
sont les proc6d£s radiographiques, et les autres des 
procedes §missifs, r6serv6s plutot aux §tres vivants et 
dans lesquels le rayonnement est produit par un corps 
emissif qu'ils ont absorb^ auparavant. Quoique les deux 

25 families de procedes soient compl^tement diff6rentes, 
la plupart des procedes de reconstruction d' images 
s'appliquent aux deux, et c'est le cas aussi de 
1' invention. 

Un rayonnement divergeant est souvent pr6f6r6 
30 pour englober plus facilement l'objet 6tudi6. Dans le 
cas d'un rayonnement d' irradiation att6nu6 par l'objet, 
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on y parvient en disposant une source Emissive k un 
foyer ponctuel prds de l f objet et un reseau de 
d6tecteurs du cote oppose de l f objet, qui sont to us 
collimates vers la source ; dans le cas d f un 
5 rayonnement 6mis par l f objet, un r6seau similaire de 
detecteurs est mis en collimation vers un foyer qui 
correspond geom^triquement & la source ponctuelle 
pr6cedente, mais a aucun objet materiel. Dans tous les 
cas f le foyer et le r6seau de detecteurs sont d§plac6s 
10 autour de I 1 objet en en prenant des vues successives, 
et pour chacune des vues, les detecteurs mesurent des 
sommes de la propriety d f attenuation ou d f Emission le 
long de lignes de collimation, qu'on appelle des 
projections de l 1 image de I 1 objet- Quand un nombre 
15 suffisant de vues a ete pris, on dispose d'un grand 
nombre de projections s 1 entrecroisant 4 travers 
l f objet. Le problSme math&matique dit d f inversion du 
syst£me permet d'obtenir la propriete d f attenuation ou 
d 1 emission, qui sert a former l f image, £ chacun des 
20 points de l 1 objet a partir des sommes de proprifetes sur 
les projections. II peut Stre repr6sent6 par un systdme 
lineaire d 1 Equations oti les valeurs connues au d6part 
sont les mesures prises par les detecteurs aux 
different s emplacements et les inconnues sont les 
25 valeurs de la propriete aux diff^rents points de 
l 1 objet. Certaines methodes de reconstruction dites 
algebriques reviennent ef f ectivement k une inversion de 
ce syst^me ; il existe cependant d'autres methodes 
dites analytiques, oix la valeur de la propri6te £ 
30 chaque point est calculee directement a partir d f une 
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combinaison math£matique des projections. L' invention 
est l'une d'entre elles. 

Un domaine privil6gi6 de la tomographie est 
l'6tude des etre vivants pour la medecine . Cornme 

5 plusieurs vues successives doivent §tre prises, l*objet 
etudi6 est susceptible de se d6placer entre elles selon 
une loi generale d 1 evolution- Les differentes vues 
repr^sentent done des etats different s de l'objet, et 
des artefacts de reconstruction ne peuvent 

10 qu f apparaltre sur l f image reconstruite. 

Des precedes pour echapper a ce probleme 
permanent consistent & prendre toutes les vues 4 la 
fois avec autant de sources et de reseaux de 
detecteurs, ou k ne prendre des vues de l'objet qu f a 

15 des 6tats identiques de celui-ci, ce qui est possible 
si son mouvement est p6riodique come le sont certaines 
activit6s physiologiques corame les battements du coeur 
ou la respiration ; le premier de ces proc§d6s est 
cependant dispendieux et le second est delicat a mettre 

20 en ceuvre convenablement . 

C f est pourquoi on a aussi entrepris 
d'incorporer des modeles num6riques du mouvement de 
l'objet refl^tant la loi g£n6rale d 1 evolution de 
l'objet dans le proc^de de reconstruction. Le module 

25 revient g6n6ralement £ une collection de champs de 
deplacements des points de l'objet aux instants 
respectifs oft les vues sont prises. II peut etre etabli 
par d f autres mesures ou par certaines hypotheses. 
Plusieurs genres de modules ont dej£ 6t6 proposes , mais 

30 ils pr6sentent en general des insuf f isances . C f est 
ainsi que le brevet US A 5 287 276, appliqu6 k l'6tude 
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de la cage thoracique, considdre les translations et 
les dilatations de son contenu, mais ce modele est 
insuffisant pour d'autres mouvements comme celui du 
cceur, qui comprend une torsion. De plus, les jeux de 

5 projections utilises pour appliquer la formule 
d 1 inversion ne sont pas connus rigoureusement, ce qui 
peut permettre de nouveaux artefacts . 

Le procede de 1 T invention pr^sente les 
avantages d'utiliser un modele plus felabore mais aussi 

10 trds simple du mouvement de l f objet, comprenant de 
nouvelles classes de mouvements, et de recourir a une. 
formule d 1 inversion nouvelle et parfaitement exacte. 

De plus, il est facile d T obtenir les jeux de 
projections necessaires et seulement n^cessaires a la 

15 reconstruction des points de l 1 image. 

Des generalisations k des situations plus 
complexes sont aussi propos6es. 

L T invention concerne sous sa forme gen£rale un 
proc6d6 de reconstruction d'une image tomographique 

20 d'un objet notamment mobile et deformable, l f image 
etant un ensemble de valeurs d'une propriety prises par 
l f objet a une collection de points de l f objet, 
comprenant l'emploi : d'un rayonnement divergent & 
partir d'un foyer et traversant I 1 objet, le foyer etant 

25 mobile autour de l f objet ; d f un module analytique de 
mobilite et de deformation de l 1 objet defini pour 
chaque position du foyer ; et d f un proc6d6 de calcul 
analytique pour obtenir lesdites valeurs k partir de 
sommes des valeurs de la propriete le long de lignes de 

30 projection menant au foyer et passant respectivement 
par les points ; caract6ris6 en ce que le moddle est 
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affine, et est une combinaison variable en cours 
d 1 acquisition , cette combinaison comprenant des 
translations , des rotations et des homotheties de 
l f objet a partir d'une origine, et en ce que le procede 
5 de calcul analytique comprend les etapes suivantes : 

- pond6ration des mesures, cette pond6ration 
etant dependante du modele analytique de 
mobilite et de deformation de l f objet ; 

- derivation des mesures ponderees suivant la 
10 trajectoire du foyer en considerant une 

direction adapt6e au module, cette direction 
etant maintenue constante, et obtention de 
mesures modifiees ; 

- retroprojection des mesures modifiees. 

15 Avant la ponderation, un filtrage des mesures 

acquises par un filtre de Hilbert est souvent applique, 
ainsi qu'il est usuel dans d'autres proced^s. La 
ponderation et la derivation s'appliquent alors aux 
mesures filtr£es. 

20 Les figures suivantes sont introduites pour 

l 1 explication de 1 T invention : 

- la figure 1 illustre la prise des mesures de 
1 1 ob j et etudie, 

- la figure 2 illustre une representation 
25 equivalente, 1 1 objet etant maintenu k un etat de 

reference , 

- les figures 3a, 3b, et 3c illustrent un cas 
complexe de deformation de l f objet f 

- la figure 4 est analogue k la figure 2 dans 
30 le cas complexe , 
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- la figure 5 illustre une prise de mesures 
avec un rayonnement conique, 

- et les figures 6, 7 et 8 sont des 
organigrammes de trois modes de realisation de 

5 l 1 invention. 

On aborde maintenant les figures . La premiere 
d'entre elles illustre la situation classique en 
tomographie d'une source S de rayonnement et d f un 
r^seau de d^tecteurs D mobile de concert sur des 

10 trajectoires circulaire s concentriques a des positions 
opposees autour d f un objet E & 6tudier. Des rayons R 
joignent la source S aux detecteurs respectifs du 
r6seau D. On consid&re ici un probldme g6om6trique 
plan f qui correspond aux conditions de mesures dites en 

15 eventail parallele oti l f etude de 1' objet E est faite 
par tranches superpos6e s. Les procedes impliquant un 
rayonnement conique sont aussi usuels et seront 
examines plus loin. L T invention a 6t6 congue pour des 
rayonneraents divergents en g6n6ral f qui englobent aussi 

20 ces derniers. 

La trajectoire T de la source S et l f ouverture 
du faisceau de rayonnement sont generalement choisies 
pour que celui-ci englobe toute la section de l f objet 
E. 

25 Consid6rons un point P de l f objet E. Un rayon 

unique provenant de la source S le traverse quand une 
vue est prise f et d'autres rayons R2, R3, etc. le 
traversent de m&ne pour d'autres positions S2, S3, etc. 
de la source S quand d f autres vues sont prises. Ces 

30 rayons passant par le point P sont utilises pour 
determiner l 1 image du point P dans le proc6de 
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d f inversion. II est utile que les rayons consid6r£s 
soieht orient6s dans des directions aussi variees que 
possible f de fagon & decrire toutes les directions de 
droite possibles pour chaque point de l T objet E. 
5 Consid6rons maintenant que l f objet E se deplace 

et se deforme. Le point P se deplace & l f int£rieur de 
la trajectoire T et les rayons R2, R3, etc. a employer 
passent par des positions respectives P2, P3, etc. 
distinctes du point P d'origine. Le probl£me de 
10 reconstruction & r6soudre est identique a la 
configuration fictive de la figure 2, ou le point P a 
et6 pris corame reference immobile et oii tout se passe 
comme si la source S suivait une trajectoire T 1 de 
forme irreguli£re et k la definition numerique notee 

15 par OS = T x (a(X)) (et a(k) dans la configuration r6elle de 
la figure 1) . On verra que cette representation fictive 
du probldme aide cl sa comprehension. D'apres 
l 1 invention, le modele de deplacement et de deformation 
qu f on choisit pour l f objet E est affine, compose par 

20 exemple d f une combinaison variable en cours 
d f acquisition de translations, de rotations et 
d'homotheties, conformement d. la formule (1) 



x 0 - r L W - A,x + B X - ^ x + [^ W J (1) 

25 oil x 0 est la position vectorielle du point P par 

rapport a une reference telle que le point O a 
l f instant de reference choisi pour mener les calculs et 

la reconstruction de l f image, x est la position du 
point P 4 un autre instant et notamment un instant de 
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prise de vue, et les coefficients a et b dependent du 
temps. X est un parani£tre gen6ral du procede, qui varie 
entre une valeur minimale X min et une valeur maximal e 
A, max ; le temps est repere par une fonction monotone 
5 de X (6ventuellement par X lui-meme) . L' instant de 
r6f£rence oii l 1 image est reconstruite correspond a X=0 . 
La propriete d' attenuation ou d' emission a calculer est 
not6e f, et pour le point P qu'on a rep6re par le 
vecteur x on la note f x (x) a un instant X quelconque et 

10 ^(xo) £ l f instant de reference, oCl x est not6 x 0 . 

La valeur d'une projection de la propriety 
mesur6e sur tout un rayon R allant de la sources a un 
d6tecteur du reseau D k travers I'objet E est donn6e 

par la formule (2) , ou a est un vecteur unitaire, et 
15 on note cette projection en fonction de ce vecteur 

unitaire a et du parametre X : 

(M* ) = £ d tf x (a(A) + ta ) ( 2 ) 

ou t est un parametre exprimant la progression 
20 sur le rayon de la projection et prenant les valeur s 
comprises dans l 1 ensemble des r6els positifs R. 

II est habituel d'appliquer un filtre de 
Hilbert aux projections dans les m^thodes analytiques, 
ce qu'on fait ici aussi, et les projections filtr^es 
25 sont notees gHA selon la formule (3) 



gnx (X, n ) — f da h H (n.a)g x (A, a) , (3 ) 
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Ou n est un vecteur quelconque, S 1 est la 
sphere - unite en bidimensionnel, Iih est donne par la 
formule (4) 

5 h H <s) =- J da isign(a)e 2i *° s . (4 ) 

— oo 

L'homme du m6tier pourra utiliser pour r6aliser 
ces formules soit une operation classique de 
convolation discrete qui fera intervenir une version 
discretisee du filtre de Hilbert, soit une operation de 
10 multiplication dans le domaine de Fourier, qui fera 
intervenir une version apodisee de la Transformee de 
Fourier du filtre de Hilbert. 

Les inventeurs ont etabli en outre la relation 
suivante (5) 

15 g^(KA T x n) = j_^P off (^ , (5) 

ou le produit matrice-vecteur A^n est la 
direction orthogonale & la droite de projection & 
l f instant X passant par le point mobile correspondant 
au point P, n est encore un vecteur quelconque et ou 
20 P oH d^signe des projections parall&les T OH (n,s), 
obtenues par rearrangement des projections 
eff ectivement mesur6es sur l'objet & l'6tat de 
reference et filtrees par le filtre de Hilbert d'apres 
la formule suivante (6) 

25 

PoH(n,s)=f ds , P 0 (n,s , )h H (s-s , ) r 

J R 

s 6tant un parametre donnant la distance £ 
l T origine. 
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La formule d T inversion pratiqu6e est alors la 
formule suivante (7) 

oft H* d6signe un vecteur unitaire orthogonal a 
un rayon reliant le point P consid6r6 eL la source S 
dans la representation de la figure 2, et ou qn. est 
donne dans la formule suivante (8) 



gFx(p,A|n)=A| det A p |g H3 (p,A|n) }=A 



10 L'homme du metier pourra utiliser une methode 

de difference finie sur deux ou trois points pour 
realiser cette 6tape de derivation. 

Cette derivation s'effectue suivant le 
parametre X de la trajectoire du foyer (la source S) . 

15 Elle est sp6cif iquement adapt6e au modele affine de 
deformation et de mobility car elle s f applique sur une 
direction a£H en consid6rant que la direction n est 
maintenue constante. II s f agit done de maintenir 
constant , non pas la direction orthogonale a un rayon 

20 acquis , mais la direction orthogonale au rayon 
equivalent de la geometrie fictive. 

A^n est la direction orthogonale 4 la droite de 
projection sl l 1 instant X passant par le point mobile 
correspondant au point P. 

25 Cette formule d' inversion tient compte des 

deformations de l'objet E et comporte, par rapport & 
d f autres formules, etablies dans des cas ordinaires, 
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une ponddration en fonction de la deformation de 
l'objet (par le determinant de A x ) et de la position de 

la trajectoire (par la norme entre x 0 et rJa(X)). Ces 
conditions peu contraignantes s'expliquent parce que, 
avec un espace affine de trans formation , les droites de 
projection de la figure 1 restent des droites dans la 
geom6trie fictive de la figure 2, de sorte que le 
probleme nurnerique peut etre resolu analytiquement . 
Dans la formule (7) 

/.&>-/ #».)*j£=£&pn— 

les homes de l f int6grale A r (x 0 ) designent un ensemble 
minimal de positions X de la source S tel que sur 
l'objet E a I'etat de r6f6rence les directions des 
droites reliant x 0 a r x (a(X)) couvrent tout l'intervalle 

15 d'un demi-tour de trajectoire sans redondance. Si on 
prend soin de limiter ainsi l f int<&grale, la formule de 
reconstruction est parfaite. L 1 integral e de la formule 
(7) pourra §tre classiquement discr6tis6e par l'homme 
de metier par une formule des trapezes. 

20 11 ©st manifeste que les deformations affines 

ne suffisent pas k mod61iser tous les objets qu f on peut 
devoir 6tudier dans la pratique. Refaisant le 
raisonnement des figures 1 et 2, on donne les figures 
3a, 3b, 3c, et 4. L'objet E 6volue de fac?on complexe 

25 entre les positions indiquees successivement aux 6tats 
Xl, \2 et A3 oil la position de la source S a aussi 6t6 
indiquee de m§me que le rayon Rl, R2 ou R3 menant au 
point P. 
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En remettant I'objet E dans un 6tat de 
reference, la figure 4 montre que les rayons notes R'l, 
R'2 et R f 3 menant au point P et correspondant aux 
rayons Rl, R2 et R3 (passant par les merries points de 
5 l f objet E) ne sont plus rectilignes . La formule (7) 

fo(*o) = 1^^ 1 h . A& Fxfr'K"*) n,est al ° rS P 1US 

valable directement. On s f y ramdne toutefois par des 
approximations/ en supposant une deformation affine. de 
10 I'objet E particuli^re k chaque point P et valable 
autour de lui. On calcule alors la matrice de 
deformation correspondant & A x de la formule (1) 



15 la formule (9) 



J-r xl O,r- 1 ^ 0 ^r xl (p,r- 1 (s 0 ) ) 



oil r xl (5i, x) et r x2 (*,, x) sont les composantes 
suivant les directions xl et x 2 principales de r(x, x) . 
Bjl n'a pas besoin d'etre calculi. Les deplacements 
20 estimes de tous les points de l f objet E sont done les 
coefficients r(k, x) qui peuvent concretement etre 
introduits dans des cartes de d6placement qui sont lues 
et exploitees au moment de 1' inversion. 

A chaque point P, il faudra encore determiner 

25 1' ensemble de projections A r (xo) qui seront utilisees 



WO 2005/001775 ^^CT/FR2004/05l»295 



13 



pour 1 T inversion, selon le principe pr6c6dent qu T un 
intervalle angulaire d f un demi-tour doit etre couvert 
par ces projections. Corame les rayons de la geometrie 
fictive R f l, R f 2 et R f 3, etc. ne sont plus rectilignes f 
5 on consid&re plutdt les directions de leur tangente au 
point P d' intersection. 

La formule d T inversion peut alors s'ecrire 
selon la formule (10) 



10 



exprim§ par la formule (11) 



det An 



fell 



p(n(r-HxoJl X)) 



15 De plus, n^fejt h) est la direction orthogonale 

a la droite acquise a 1' instant X passant par r -1 (x 0 ) ; 
pp(r _1 fc 0 )l A.)) est un facteur lie a la deformation des 
longueurs sur les rayons, d'apres la formule (12) 

20 Enfin, l(xo - I\(a(A,)) est la distance du point P £ 

l 1 image de la source sur l'objet de r6f6rence, le long 
du rayon virtuel R f . 

Au detriment de la quality, mais au profil 
d'une plus grande simplicity, 

25 - le terme L(x 0 -r(a(l))) peut etre remplace par 
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-r,(a(x)J 

detAjr-^xo)) 



- le terme 



par 1 . 



Le module affine peut §tre defini d f apr6s les 
deformation r6elles de l'objet par une approximation 
5 d'apr&s un critere d 1 approximation tel le critere des 
moindres carr6s, le critdre de la minimisation de la 
norme L 1 et L 2 , eventuellement completes par une 
regularisation sur le gradient ou le laplacien, 

On s'est int6ress6 jusqu'ici a des 

10 reconstructions de l T objet E par tranches, dans des 
conditions en eventail parallele. Le proced6 pr6c6dent 
peut etre 6tendu avec des acquisitions par un 
rayonnement conique dans les conditions classiques de 
la figure 5, oil l'objet E est repr6sent6 dans son 

15 enti6ret6 tridimensionnelle, de meme que le reseau de 
d6tecteurs D, qui est compose d'une s6rie de couches de 
detecteurs semblables chacune £ celle des figures 
pr6c6dentes ; un rayon R peut §tre exprim6 par trois 
param&tres, 4 savoir la position X de la source S, un 

20 angle ■ y que fait le rayon R dans le plan de la 
trajectoire T par rapport & l'axe central X du 
faisceau, et une cote q reperant la couche de 
detecteurs k laquelle le rayon aboutit. L f angle y 

remplace le vecteur de direction a precedemment 
25 utilis6 dans les formules par commodite- 
La formule suivante (13) 
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donne 1' inversion qu'on entreprend alors dans ces 
conditions g6om6triques particulieres et en faisant une 
approximation appelee de Feldkamp courante dans l f art 
5 pour traiter les projections coniques. Les mesures sont 
ici multiplies par cos A, qui est la pond6ration 
class ique compensant le caract^re oblique du rayon 
acquis. Cet angle A est donne sur la figure 5 et 
repr6sente 1 T angle du rayon R avec le plan de la 

10 trajectoire T. De plus, I*(xo - r x (a(X)) est la distance de 

x k l 1 image de la source S sur I'objet E & l f 6tat de 
r6f6rence le long du rayon virtuel R F / en projection 

sur le plan de la trajectoire T, et ^(r-^xoU) est la 

direction orthogonale £ la droite acquise & 1" instant X 

15 passant par T^xo) sur ce plan. Enfin, on a vu que 

q(k, r^fcO) est la cote sur le reseau de detecteurs D de 

la droite, acquise a l 1 instant X, passant par le point 

mobile correspondant a xo . 

Un proc6de particulier, d6crit dans le brevet 

20 frangais d f abord d6pos6 sous le numero 01 07918, peut 
encore §tre appliqu6 ici moyennant quelques 
transformations. Ce proc6d6 s f applique a des 
acquisitions rapides nien§es autour de l f objet en 
faisant plusieurs tours et consiste en une 

25 reconstruction de I'objet anim6 en evaluant son contenu 
par des blocs de projections prises sur une fraction de 
tour seulement de la source. Les inversions de 
projections menees sur les blocs produisent des sous- 
images qui sont incorrectes puisqu'elles ne comprennent 

30 qu'une partie des mesures, mais sont obtenues trds 
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rapidement et permettent d'obtenir avec une m§me 
rapidite la loi de d6placement ou de deformation de 
1" objet en comparant des sous images homologues, prises 
pour des memes positions de sources espacees d f un tour 
5 complet ou d T un demi-tour ; les sous-images sont 
reconstruites & un instant de reference pour chacun des 
groupes de blocs et finalement combinees entre el les 
pour donner 1 1 image complete de l f objet. Ici, les 
limites des blocs sont donnees dans la g6om6trie 

10 virtuelle. Le proc6d6 n'est pas par ailleurs modifie. 

Enfin, il est possible d T am§liorer le procede 
pour l T 6tude d f un objet soumis a des phenomdnes 
p6riodiques. On s^lectionne alors comme blocs d f etude 
seulement ceux qui sont pris a des phases du ph^nomene 

15 analogues a celles de l 1 instant de reference. 

On pourra se reporter aux organigrammes des 
figures 6, 7 et 8 pour finir de prendre connaissance de 
l f invention. lis d6taillent trois modes de realisation 
de 1* invention, oH la reconstruction est faite 

20 respectivement avec compensation du mouvement ; 
compensation du mouvement et compensation temporelle ; 
et compensation du mouvement, compensation temporelle 
et prise en compte d'une periodicity de 1' evolution de 
l f objet. Ce qui a 6te decrit auparavant est applique 

25 dans ces proc6des. Les voxels sont 6videmment les 
points de l f image de 1 T objet qu'on considere dans 
l 1 image reconstruite. 
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REVENDICATIONS 

1. Proced6 de reconstruction d'une image 
tomographique d'un objet notamment mobile et 
5 deformable, 1' image 6tant un ensemble de valeurs d'une 
propri6t§ prises par l 1 objet & une collection de points 
de l f objet, comprenant l'emploi : d'un rayonnement 
divergent a partir d f un foyer et traversant l f objet, le 
foyer etant mobile autour de 1 T objet ; d'un module 

10 analytique de mobilite et de deformation de 1' objet 
defini pour chaque position du foyer ; et d'un proc<§d6 
de calcul analytique pour obtenir lesdites valeurs a 
partir de sommes des valeurs de la propri6t6 le long de 
lignes de projection menant au foyer et passant 

15 respect ivement par les points ; caracterise en ce que 
le modele est affine r et est une combinaison variable 
en conrs d' acquisition, cette combinaison comprenant 
des translations , des rotations et des homoth^ties de 
1' objet a partir d'une origine, et en ce que le procede 

20 de calcul analytique comprend les etapes suivantes : 

- pond6ration des mesures, cette pond6ration 
etant dependante du module analytique de 
mobility et de deformation de 1' objet ; 

- derivation des mesures pond6r6es suivant la 
25 trajectoire du foyer en considerant une 

direction adapt6e au modele, cette direction 
6tant maintenue constante, et obtention de 
mesures modifiees ; 

- r6troprojection des mesures modifiees. 
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2. Proc6d6 de reconstruction d'une image 
tomographique selon la revendication 1, caracterise en 
ce qu'il comprend une etape de definition d f intervalle 

d 1 integration donnant fo(xo) de fagon que les lignes de 
5 projection aient des directions s'6tendant sur un 
intervalle de N demi -tours consid^re dans un 6tat de 
I'objet oil l 1 image est reconstruite, N 6tant sup6rieur 
ou 6gal a 1. 

10 3. Procede de reconstruction d T une image selon 

l'une quelconque des revendications 1 ou 2, caracterise 
en ce qu f il est applique avec un rnodele analytique de 
mobility et de deformation particulier k chaque point 
de l'objet. 

15 

4. Procede de reconstruction d f une image selon 
la revendication 3, caracterise en ce que la partie 
lineaire affine du modele analytique particulier 
correspond & 1 T approximation locale affine d'une 

20 famille de trajectoires passant par des points situ6s 
au voisinage de chaque point, 

5. Procede de reconstruction d'une image selon 
les revendications 2 et 3, caract6ris6 en ce que les 

25 directions des lignes de projection sont considerees au 
point de l f objet oil le proc§d§ est appliquds. 

6. Proc6d6 de reconstruction d'une image selon 
l'une quelconque des revendications 1 & 5/ caracterise 

30 en ce que la reconstruction s'effectue d f abord par 
sous-images construites avec des lignes de projections 
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ayant des directions comprises dans les intervalles 
angulaires plus petits qu ? un demi-tour, en ce qu'il 
coraprend une 6tape de determination des intervalles 
angulaires auxquels appartiennent les lignes de 
5 projections et en ce qu'il comprend une etape de 
combinaison des sous -images apres correction de la loi 
devolution g6n6rale de l'objet entre les positions 
associees a differents instants de r6f6rence associes cL 
chaque intervalle angulaire. 

10 

7. Procede de reconstruction d'une image selon 
l'une quelconque des revendications prec6dentes, 
caract6ris6 en ce que la loi d' evolution g6n6rale de 
l f objet est periodique et en ce que les instants de 

15 reference sont choisis pour la meme phase des periodes 
de mouvement de l f objet. 

8. Proc6d6 de reconstruction d f une image selon 
l f une quelconque des revendications pr6cedentes, oil le 

20 rayonnement est conique, caracterise en ce qu T avant 
l f 6tape de filtrage des mesures acquises, le proc6d6 
comprend une etape de ponderation des mesures acquises 
adaptee a une acquisition des mesures en g6ometrie 
conique et en ce que la retropro j ection est effectu6e 

25 selon la g6ometrie conique en tenant compte du modele 
analytique de mobilite et de deformation de l'objet. 

9. Proced6 de reconstruction d f une image selon 
I'une quelconque des revendications prec6dentes, 

30 caract6ris6 en ce que la ponderation des mesures et 
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pr6c£d6e par un filtrage des mesures acquises par un 
filtre de Hilbert. 

10. Procede de reconstruction d'une image selon 
5 l'une quelconque des revendications pr6c§dentes, 
caract6ris6 en ce que le modele affine est obtenu par 
approximation de deformations reelles de l T objet selon 
un critere d f approximation. 



WO 2005/001775 




T/FR2004/050295 



3/6 





FIG. 5 
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